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Abstract 



We give an 0(n 1,5 logn) time algorithm for finding the maximum flow in a directed planar graph 
with multiple sources and a single sink. The techniques generalize to a subquadratic time algorithm for 
bounded genus graphs. 

1 Introduction 

In general graphs, multiple source flow problems are reduced to single-source problems but connecting 
a super source to the sources with infinite capacity arcs. In planar graphs this reduction destroys the 
planarity. Using the maximum flow algorithms for general graphs, a multiple-source, single (or multiple) 
sink max flow in directed planar graphs can be solved in 0(n 2 log n) time using Goldberg and Tarjan's 
preflow push algorithm [7] or 0(n log n log U) where U is the maximum edge capacity using Goldberg 
and Rao's binary blocking flow algorithm [5]. In this paper we give an 0(n 15 log n) algorithm for the 
problem in planar graphs by combining preflow push and augmenting path algorithms. 

In planar graphs, multiple source and sink flow problems have been studied by Miller and Naor, 
giving subquadratic-time algorithms for the case when the sources and sinks are on a common face and 
for the feasibility problem (the amount of flow out of every source and into every sink is known) [TT] . The 
maximum single source, single sink flow in a directed planar graph can be found in 0(n log n) time pQ. 



We are given a directed planar graph G with arc capacities a set of source vertices S and sink vertices 
T. A flow is an assignment of values to arcs not exceeding the capacity such that the flow entering 
a non-source, non-sink vertex is equal to the flow leaving the same vertex. A flow is maximum if it 
maximizes the amount of flow leaving S and (equivalently) there is no residual path from any vertex in 
S to any vertex in T. An arc is residual if the flow is less than the capacity. The reverse of an arc is 
residual if there is flow on the arc. For more formal definitions, see pQ. 

A preflow is a flow that allows excess inflow at vertices that are not sink vertices. A maximum preflow 
is a preflow that maximizes the flow into the sinks. We will use the following lemma, which holds for 
general graphs. The forward direction follows from the definition of maximum preflows. The reverse 
direction follows from the Max Flow, Min Cut Theorem: if there are no residual paths from sources or 
vertices with excess inflow to sinks, then there is a saturated cut separating the sources from the sinks 
and the preflow cannot be increased (see [6]). 

Lemma 2.1. A preflow is maximum if and only if there is no residual path from a source to a sink or 
from a vertex with excess inflow to a sink. 

We also assume for this draft that the reader is familiar with recursive subdivisions of planar graphs. 
A piece is a subgraph resulting from this decomposition. See [3] for details. 



2 Definitions 



i 



3 Algorithm for a piece not containing sinks 



Our multiple-source, single sink algorithm (Section |4| is a recursive algorithm using a recursive decom- 
position based on cycle separators. In this section, we present a solution to the non-trivial subproblem in 
which we need to find a maximum flow from a set of sources in a piece (of the recursive decomposition) 
to sinks outside the piece. More formally: 

Input: a piece P with a constant number of holes, a set Sp of sources in P, and a set Tp of sinks none 
of which are internal vertices of P. 

Output: The max flow from Sp to Tp. 

The algorithm proceeds in phases, mimicking a preflow algorithm in two pushes: one to the boundary 
of the piece (Phase 1) and one from the boundary to the sinks (Phase 2). In a final phase, the (maximum) 
preflow is converted to a maximum flow. 

Phase 1: Find the max flow in P from Sp \ dP to dP. 

Phase 2: For each vertex p in dP, if p is a source, compute a max pTp-flow. Otherwise, push as much 
flow from p to Tp as possible while not exceeding the inflow to p. 

Phase 3: Convert the above maximum preflow into a maximum flow. 



3.1 Correctness 

We need the following lemma, which holds for general graphs. 

Lemma 3.1. In a graph, let A, B be a partition of the vertices such that all the arcs from A to B are 
saturated and the sinks t are vertices of B. Augmenting a preflow cannot make this cut residual. 

Proof. By assumption, there cannot be an augmenting path from a vertex in A to a sink t, so consider 
an augmenting path P from a vertex b G B to t. Assume for the sake of contradiction that P makes the 
cut residual. Then it would have to intersect A. Since b,t £ B, there is an edge (u, v) on P such that 
u G A and v G B. But by assumption, this edge is saturated, which is a contradiction. This shows the 
lemma. □ 

We show that after Phase 2, a maximum preflow has been computed. Let S' P = Sp \ dP and t G Tp. 

After Phase 1, there are no residual s-Xo-t paths for any s G S' P : such a path would have, as a prefix, 
a residual s-to-p path for a vertex p G dP; this would contradict the maximality of the S'p-to-dP flow 
At the end of Phase 1, there is a saturated Sp — dP cut Cg' p (separating S' P from dP). By Lemma 
the augmentations in Phase 2 cannot introduce a residual s-to-t path for any s G S' P . 

It remains to show that at the end of Phase 2 there are no residual paths to t from a boundary vertex 
p which is either a source or which has excess inflow; by Lemma |2.1| the flow at the end of Phase 2 is a 
maximum SpTp-preflow. If p is a source, we saturate all edges in a pt-cut when finding a max pTp-flow 
in Phase 2. By Lemma |3.1| this cut will stay saturated. Now, assume that p has excess inflow f v at the 
start of Phase 2. There are two cases. If we are able to route f v units of flow to t, then p no longer 
has excess inflow. If we are unable to route f p units of flow to t, there is a saturated pt cut C p that is 
a witness to there being no residual p-Xo-t paths. Lemma |3.1| implies that augmenting flows from other 
boundary vertices cannot make this cut residual. 
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3.2 Running time 

We analyze the running time of Phases 2 and 3. Phase 1 will be computed recursively; the analysis of 
the recursive algorithm is in Section [4] 

Phase 2 can be solved in 0(|9P|nlogn) = 0( \J\P\n log n) time by computing a maximum single- 
source, single-sink flow (in the entire graph) for each p G dP and t G Tp. If p is a source then compute 
the max pt-flow algorithm for each t G Tp. If p is a non-source vertex with excess / from Phase 1, 
augment the graph with a source s connected to p by an arc with capacity / and compute the maximum 
st flow for each t G Tp. Each max flow computation takes 0(n log n) time in planar graphs. Hence, 
Phase 2 runs in 0(^/\P~\n log n) time. 
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Maximum preflows can be converted into maximum flows in sparse graphs in 0(n log n) time (see 
Goldberg and Tarjan [HE]; the conversion is only made explicit in the earlier conference version): first 
eliminate cycles of flow (0(n log n) time) and then eliminate excess inflow from vertices by processing 
them in reverse topological order in the flow's support graph (O(n) time). In fact, our maximum preflow 
is already acyclic, since the flows computed in Phase 2 are acyclic; converting to a max flow will, in fact, 
only take linear time. 

The total time required for Phases 2 and 3 is therefore 0(\/\P\n log n). 



4 Multiple-source, single-sink algorithm 

In this section, the overall max flow algorithm is described. Since we apply recursion, we need to 
consider a slightly more general problem than multiple source, single sink max flow in order to ensure 
that subproblems are of the same form: we need to find a max flow from a set of sources to at most a 
constant t number of sinks. 

Let us regard the entire graph as a piece P with no boundary vertices. We define each sink to be a 
boundary vertex of P. Each such sink now defines a degenerate hole of P and P contains no sinks in its 
interior. By assumption, the number of sinks is bounded by t so the number of holes of P is at most t. 

Next, we obtain a subdivision of P into a constant number p of subpieces; each subpiece has at most 
c p n vertices and y/c p n boundary vertices (constant c p < 1). This is done with the recursive r-division 
algorithm of Frederickson (see Lemmas 1 and 2 in [3]) but with Miller's cycle separator theorem [10] 
instead of the separator theorem of Lipton and Tarjan [9]. Through the recursion, we ensure that the 
number of holes in each piece is bounded by t— 1. If we pick t to be a sufficiently large constant, we can 
ensure this bound with the approach introduced by Fakcharoenphol and Rao 

We now run the algorithm of the previous section on each subpiece P' . To solve Phase 1 for P' , we 
define a new planar graph P" from P' by introducing a super sink for each hole and adding an edge of 
infinite capacity from each boundary vertex of that hole to the super sink. The same is done for the 
external face of P' if it contains boundary vertices. Now, solving Phase 1 for P' corresponds to finding 
a max flow from the sources in P" to its super sinks. Since P' has at most t — 1 holes, P" has at most t 
super sinks and so we can recurse on P" to solve this problem (note that only the super sinks of P" and 
not the sinks of P belonging to P" are regarded as sinks in the recursive call). 



4.1 Running Time 

Phases 2 and 3 run in 0{n 3 ^ 2 logn) time. 
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As described in Section 

We define a recurrence relationship to bound the time for Phase 1 and the entire algorithm. For n 
larger than some constant, the size of each piece in the division (including the super sinks added to the 
piece) is at most c' p n, c p < c' p < 1. 

We repeat the recursion until each piece has size at most some constant r. Let N be the total size 
of all pieces at the leaves of the recursion tree. Using ideas of Frederickson [4] (see Lemmas 1 and 2 and 
their proofs in that paper for details), N = kn for some constant k. 

In the analysis, we assume that the total size of pieces at any recursion level is N by regarding the 
vertices in pieces at the leaves of the recursion tree to be present in every level. This allows us to regard 
the pieces in any level as being pairwise vertex-disjoint by counting a vertex according to its multiplicity, 
once for each piece to which it belongs. Furthermore, we may assume that no new vertices (i.e. super 
sinks) are introduced in the recursive steps since all N vertices (including all added super sinks) are 
present in any recursion level. Finally, we may assume that the size of a piece in a division is exactly 
c p n since larger pieces will only increase running time. 

From these assumptions, it follows that we may restrict our attention to the case p = l/c' p and the 
recurrence relation for the running time of our algorithm becomes 

T(N) < pT(N/p) + c'N 3/2 log AT, 

where c' is a constant and the second term is the total time for Phases 2 and 3 and for computing the 
division of the graph. We will show that T(N) < cN 3 ^ 2 log N for some constant c (which we may assume 
is true for small N). For the induction step, we have 

T(N) < pc(N/p) 3/2 log N + c'N 3/2 log N = (c/^p + c')N 3/2 log TV. 
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The desired time bound is achieved by setting c= c'/(l — 1/^/p) > 0. 



4.2 Correctness 

Correctness of the overall algorithm follows from Theorem |4.1| that is, the algorithm is correct over all 
pieces because we can consider the sources in any order, and we always fully saturate the flow from the 
source we consider. 

Theorem 4.1. Given an instance of a multiple source, multiple sink maximum flow problem, one can 
find the maximum flow by: 

iterating over the sources s in any order 
iterating over the sinks t in any order 

saturating a maximum st-flow in the residual graph 

Proof. We prove by induction over the inner and outer loops of the algorithm. Let si,ss, ... be the 
arbitrary order that the sources are iterated over. 

Outer loops: Say we are currently considering source Sj. Assume that there are no residual paths 
from source Sj for any i < j to any sink. We will show that no residual Sj-to-t path can, by way of 
augmentation, introduce an Si-to-t' residual path. By the inductive hypothesis, the Max Flow, Min 
Cut Theorem, and submodularity of cuts, there are non-residual non-crossing cuts Si, Si and S' t , S'i with 
Si S Si n S'i and t £ Sj PI S'i. Let R be a residual Sj-to-t path. If R could introduce a residual Si-to-t' 
path, then it must be that Sj G Si and t S Si. By the submodularity of cuts, it must be that Si C Si. It 
follows that Si, Si is a non-residual cut that also separates Si and t' and augmenting R cannot introduce 
a residual Si-to-t' path. 

Inner loops: Suppose we are currently considering source sg. Let ti, ta> • - • be the arbitrary order that 
the sinks are iterated over for source St- Suppose we are saturating the max sitj flow. By the above 
argument, doing so will not introduce a residual source-to-sink path for any source Sk, k < £. Assume 
that there are no residual se-to-ti paths for any i < j. By this inductive hypothesis and the Max Flow, 
Min Cut Theorem, there are non-residual cuts Si, Si such that si £ Si and ti € Si. If there is a non-zero 
se-to-tj flow in the residual graph, then it must be that tj £ Si for every i < j. It follows that no se-to-tj 
flow can leave Si. □ 

Note that saturating an arbitrary order of the source, sink pairs will in general not result in a maximum 
flow, see Figure [l] 




Figure 1: A counterexample to arbitrarily saturating source sink pairs: the maximum {s, s'}-to-{£, t'} flow 
has value 3. If the source, sink pairs are maximized in the order (s, t), (s' , t'), (s, t'), (s' , t), the value of the 
flow found is only 2. 



5 Bounded Genus Graphs 

Since our algorithm does not rely too much on planarity, it is easy to generalize it to bounded genus 
graphs. For such graphs, we can use the 0(n log 2 nlog 2 C) time max st-flow algorithm due to Erickson, 
Chambers and Nayyeri [2] (C is the sum of capacities). To begin the decomposition, we may first 
planarize the graph with a set of g simple cycles, each of length 0(i/n) (e.g. construction by Hutchinson 
and Miller [8]). This will allow us to continue with a decomposition based on planar separators and 
introduce only O(g) holes in which we can embed a super sink while maintaining planarity as required 
for Phase 1. This gives an 0(n 3 ^ 2 log 2 nlog 2 C) time algorithm; of course, using the binary blocking flow 
algorithm would be asymptotically faster. 
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